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We use the Hubbard type model to describe the coexistence between superconductivity (SC) and 
ferromagnetism (F). Our Hamiltonian contains single-site and two-site interactions. All inter-site 
interactions will have included the inter-site kinetic correlation: {cf a Cj a ), within the Hartree-Fock 
approximation. To obtain the SC transition temperature Tsc and Curie temperature Tc we use 
the Green's functions method. The numerical results show that the singlet SC is eliminated by F, 
but the triplet SC is either enhanced or depleted by F, depending on the carrier concentration and 
direction of a superconducting spin pair with respect to magnetization. The kinetic correlation is 
capable of creating superconductivity. We find that the ferromagnetism created by change of the 
bandwidth can coexist with singlet superconductivity. In the case of triplet superconductivity the 
ferromagnetism creates different critical SC temperatures for the A\ and A2 phase (the pair's spin 
parallel and antiparallel to magnetization, respectively). 



PACS numbers: 71.10.Fd, 74.20.-z, 75.10.Lp 



1. Introduction 

In the last few years materials have been found where superconductivity (SC) coexists with ferromagnetism (F) 
within the same electron band. These two cooperative phenomena are mutually antagonistic because supercon- 
ductivity is associated with the pairing of electron states related to time-reversal, while in the magnetic states the 
time-reversal symmetry is lost, and therefore there is strong competition between them. However, Ginzburg [l[ has 
pointed out the possibility of this coexistence under the condition that the magnetization is smaller than the ther- 
modynamic critical field multiplied by susceptibility of a given material. Matthias and co-workers [2| demonstrated 
that a very small concentration of magnetic impurities was enough to completely destroy superconductivity when 
ferromagnetic ordering was present. 

Superconductivity and ferromagnetism was simultaneously observed in the late 1970s in the intermetallic systems, 
e.g. ErRh4B4 Q and HoMoeSg [4]. These systems have a Curie temperature Tc much lower than the critical su- 
perconducting temperature T$c- The domain of coexistence is in a very limited narrow temperature range. These 
compounds also have the second critical temperature Tsci < Tc below which the superconductivity disappears and 
they become only ferromagnetic. In both: ErRh 4 B 4 and HoMogSs, the ferromagnetism and superconductivity are 
created on different atoms, e.g. in ErRh 4 B4 ferromagnetism is carried by 4/ electrons of Er atoms and the SC by Ad 
electrons of Rh atoms. 

There is an intermetallic compound exhibiting a coexistence of SC and F, this is ErNi2B2C which is characterized 
by the long range magnetic order associated with magnetic Er ions. This compound has Tsc ~ 10-5 K; below TW ~ 6 K 
superconductivity coexists with antiferromagnetism, and below Tsc ~ 2.3 K SC coexists with ferromagnetism [6j. The 
magnetic moment of this compound is very weak and is equal to 0.35^_e/Er atom at temperature 1.8 K, which is 
about 1/22 of the saturation moment Q. 

There are also results showing the coexistence between magnetic ordering and superconductivity in a family of 
hybrid ruthenate-cuprate compounds such as RuS^GdC^Os- This compound exhibits a ferromagnetic order at 
a rather high Curie temperature Tc = 133 — 136 K, and becomes superconducting at a significantly lower critical 
temperature T sc = 15 -40K 

Recently the so-called ferromagnetic superconductors were discovered, which exhibit simultaneously the ferromag- 
netic and a spin triplet superconducting phase modified by external pressures. To these ferromagnetic superconductors 
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we include UGe2 0, j URhGe 11 1 , and ZrZn2 [l2j ■ Their Curie temperature is much higher than the critical su- 
perconducting temperature, and the superconductivity exists only in the presence of ferromagnetism. Both states: 
ferromagnetism and superconductivity are formed by the same electrons. In UGe 2 and URhGe these are 5/ electrons 
on U atoms, and in ZrZn 2 these are 4c? electrons of Zr atoms. 

The theoretical possibility of ferromagnetism coexisting with the triplet parallel spins superconductivity was sug- 
gested by Fay and Appel for ZrZn 2 [HI , while the coexistence of F with singlet superconductivity was theoretically 



developed by Fulde, Ferrell [14l |. and Larkin, Ovchinnikov [15(. Further theoretical development took place after 



finding experimental evidence for the coexistence of triplet SC with F (see e.g. HI, l3|) or singlet SC with F (see e.g. 

In this paper we propose the model describing coexistence of ferromagnetism and superconductivity which is based 
upon the kinetic interactions At and t ex , and the inter-site kinetic correlation I a = (c+ Cjo-). It was shown 2(], 21, 22| 
that the ferromagnetic order can be created in two ways: (i) due to mutual shift of majority and minority spin bands, 
as in the Stoner model, (ii) due to a change in the bandwidth of the spin bands. The first case leads to decrease of 
potential energy of the system, and the second one leads to lowering of the kinetic energy. The kinetic interactions 
give rise to both of these mechanisms in creating the ferromagnetic state. They were also used in describing the 
superconductivity [H, 0, [H, [5^, US Ell • Hirsch dSHI and Micnas [HI have shown that the interaction At leads 
to s-wave singlet superconductivity in a less then half-filled band system. Aligia and co-workers have demonstrated 
"28L 130] ] that the interaction t ex can give rise to both: singlet superconductors with s-wave and d-wave symmetry 
28| and triplet superconductors with p-wave symmetry (30j . 
The paper is organized as follows. In Section 2 we have put forward the model Hamiltonian and we analyzed this 
Hamiltonian using Green's functions technique. In Section 3 we examine the possibility of the coexistence between 
singlet superconductivity and ferromagnetism. We show the dependence of Curie temperature and superconducting 
critical temperature on carrier concentration. A similar analysis for triplet superconductivity is presented in Section 4. 
In Section 5 we present the influence of external pressure on the Curie temperature and the superconducting critical 
temperature for the triplet superconductors. Section 6 contains conclusions covering results presented in this paper. 




2. Model Hamiltonian and Green's functions technique 

The experimental results for ZrZn 2 show that the magnetic moment is 0.17 /is, which is a small value. This 
suggests that the ferromagnetic superconductor ZrZn 2 is a weak itinerant ferromagnet. It will allow a description of 
this compound by the extended Stoner model. In this model we include the on-site Coulomb repulsion and the inter- 
site interactions which modify the electron bandwidth. In the superconducting state we include the experimentally 
observed spin triplet superconductivity and also the spin singlet superconductivity. Our Hamiltonian can be written 
in the following form [3l|, H3] 



H = - [t - At (n,_ CT + nj- a ) + 2t ex ri i - (T hj- lT ] c+ a c, J(y - fi ^ h ia + — ^ n ia n t ^ a + - ^ cf. cf^Ci^Cj^, (1) 

<ij> ia ia <^j> 

aa 

where cf a (ci a ) creates (annihilates) electrons with spin a on the lattice site i, hi a = c^c^ is the electron number 
operator for electrons with spin a on the lattice site z, U is the on-site Coulomb repulsion, J is the inter-site exchange 
interaction, and ij,q is the chemical potential. Kinetic hopping and exchange hopping interactions: At and t exi 
according to [32] , can be expressed as 

At = t-t!=t(l-Sx) and t ex = — ^ - t a = t ^ -, (2) 

where Si and 5* 2 are the hopping inhibiting factors 



tx/t = S x and t 2 /h = S 2 , (3) 

t\ and t-i are the electron hopping energies in the presence of one and two electrons with opposite spin, respectively. 

After applying the modified Hartree-Fock approximation (H-F), in which the inter-site kinetic correlation {cf a Cj a ) 
is included for all inter-site interactions, we obtain the following simplified Hamiltonian (see §6 of Ref. 32]) 
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H = H o + Hosp + Hesp- (4) 
The Hamiltonian Hq is the kinetic energy with added molecular field 

Ho = - ta eff c ta c 3<y -J2(^o- M a ) h ia , (5) 

<i}> ia 
a 

where the effective hopping interaction is t G e ^ = W , the bandwidth modification factor has the form 
1 



b° = 1 

t 



2Atn_ ff -2i ex (n^-/2 (T _2/ CT /_ CT -|A | 2 -|Af| 2 + |A^| 2 ) + J (/„. + /_„)] , (6) 

and the spin-dependent modified molecular field is given by 

M a = Un- a - zJ n<y + 2zAtI a - 2zt ex (2J_ CT n CT + Af,A* + Af*A ) , (7) 

where z is the number of nearest neighbors. 

The electron occupation number used above is n a = (c^ r Ci a ') and the Fock's parameter, proportional to the kinetic 
energy, is given by I a = (c+ c Jcr ). The bandwidth change factor V includes the following averages 

A = ~ Y a (Ci-aCta), Af = i Y ° ( c *i+h-aCia) , A£ = (c t+h , a C l<7 ) , (8) 
a (7 

where h = j — i is the difference of lattice indices for the nearest neighboring atoms in a given direction, Ao and 
A^ are the single-site and inter-site singlet superconducting parameters, and A^ is the inter-site equal spin triplet 
superconducting parameter. 

The next two terms of the Hamiltonian ((4]), Hosp and Hesp, are related to the opposite spin pairing (OSP) and 
equal (parallel) spin pairing (ESP), respectively. The OSP term describes the singlet SC with the total spin being 
and the triplet SC with total spin 1 (in units of Ti) and its projection being (see (33|). This term is given by 

Hosp = J2 ai (c+c+ + /i.e.) + ^ a 2 (c+c+ - c+c+ + h.c) + ^ «3 (c£c+ + c+c+ + h.c) , (9) 

i <ij> <ij> 

where 

oi = [U + 2zt ex (I T + /;)] A + 2z (At - t ex n) Af , (10) 



a 2 = (At - t ex n) A + ^ [J + 2t ex (J T + J;)] A s h , (11) 
«3 = i[-J + 2t e ,(/ T +/ i )]A^, (12) 

and 



is the inter-site opposite spin triplet superconducting parameter. 

The ESP term describes the triplet SC with total spin 1, its projection being ±1 (33|, and is given by 
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h esp = \ £ (- J - 4 *ex/- CT ) A£ (c+c+ + h.c). (14) 

<ij> 

a 

Transforming Hamiltonian (J4j) into the momentum space we obtain 

H = J2^k, - £ [(Ali + A&) c +c+_ kl + h.c] ]T (AZc+ a c+_ ka + h.c), (15) 

ku k ka 

where 

4 = e k b a + M CT - Mo (16) 

is the spin dependent modified dispersion relation, e k is the unperturbed dispersion relation which has the following 
form 

£ fc = -<$> ik - h , (17) 

h 

where h — r, is pointing to the nearest neighbors lattice sites. 

The singlet superconductivity energy gap A^i is the Fourier transformation of the first two terms in the Hamiltonian 
([9]). It is symmetric, A_ feg = A k g, and is given by the relation 

A&^do + diAS+daAg), (18) 

where 

d = - [U + 2zt ex (J T + /;)] A - 2z (At - t ex n) Af , 

di = - (At - t ex ii) , (19) 
d 2 = -±[J + 2t ex (Ii+I l )}. 

The momentum independent gap parameter above is given by a Fourier transformation of Eq. ([H]) to momentum space 



ka 

and the momentum dependent gap parameters are given as 

A« = A J2 + ^ k - h ) » AS = £ Af (e 4k ' h + e"**) . (21) 
ft ft 

The opposite spin triplet energy gap A£j> is the Fourier transformation of the third term in the Hamiltonian (J9]), 
it is antisymmetrical, A_ feT = — A^, and is given by 

A^ = d 3 A kT , (22) 

where 

d 3 = ±[J-2t ex (I< [ + I i )], (23) 
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and 



h 



T I „ik h -ik-h> 



(24) 



The triplet superconducting ordering parameter for parallel spins is denoted by A k and it is the Fourier transfor- 
mation of A a h from Eq. ([8]). It is given by the relation 



where 



A% = d^ kT , 



di = -{J + Uexl-a) , 



(25) 



(26) 



and 



A a \ A a a ik-h 

kT=Z^ A h e ■ 



(27) 



The Hamiltonian (fT5|) will be analyzed using Green's functions technique. The Green's function ((A; B)} £ satisfies 
the following equation of motion 



e((A;B)) £ = ([A,B} + ) + (([A,H}_:B)) s . 
Using Hamiltonian (|15|) in the equation above we obtain a set of equations, which can be written as 



(28) 



T 



-2(Al fc ) 
.Tl \* /aTI 







-(A u 



A S 



Y (a!U 



2Al 



.Tl 



^kS 



A u \ 

^kT \ 



2A 



-2(Ai fc )* 



e + e 



G(k,e) 



— k 



where 1 is the unity matrix and the Green's functions matrix has the form 



G(*,e) 



/ (( c feT; c feT)>£ (( c fcT;4i))s {{ck]\c_ k] )) e ((c fcT ;c_ fei )) e \ 

(( c fei; c fet))e (( c ki;ct))e (( c ki;c-ki))e ((c fei ;c_ fei )) e 

\ i( C -kl> c kT»« ((ctfc;;Cfc|»e (( c lfei; c -fcT))e ((cl fci ;c_ fei » £ / 



Using the symmetry properties of the order parameter we can rewrite Eq. (|29[) in the following form 



(29) 



(30) 



e — ei 



2Al 

.Tl 



a;.l. -- Ali \ 



2(a1; 



V(A 



kS> 



(A 



Tl \ 
kT> 



2(Ai)* 



-A^ + A fcT 




Tl 

fcS "T" ^fcT 

2Ai 







£ + £ 



G(k,e) = i. 



(31) 



In further analysis we will consider separately, the coexistence of ferromagnetism with the singlet, and with the ESP 
triplet superconductivity. 



3. Coexistence of ferromagnetism and singlet superconductivity 
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In this case the energy gaps in Eq. ([31]) are: A£ = 0, A^jj, = 0, and ^ is given by Eq. [[IS]). This 

assumptions simplifies the Green's functions G (k, e). Using in Eq. (JSJ the Zubarev's relation [UlSHl for the average 
of the operators' product, and the Green's functions G (k, s) from Eq. (|3l"|) . we calculate the on-site superconductivity 
parameter Ao 

A = ^^2<T (Ci-trCiv) = ojy E ^ ( C -k-<rCka) = [ f ( e ) Im ((( c fcT > c -kl))e - ({Ckf, C_fc T )) e ) <fe, (32) 



and the inter-site superconductivity parameter A 



£7 k(T k 

(33) 

To solve the above equations we obtain functions ((cki; c_^)) e and ((cfcj,; c~ k ^)) £ from Eq. (|3ip . They have the 
following forms 



A 



«c fcT ;c_ fci )) E = -- — SS_ (34) 



- I A Ti 



and 



((c fei ;c_ feT )) e = - — - T . (35) 

(-4) (W^)-(O 

Inserting these functions to Eqs. (|32)) and (|33p we obtain 

A o = -7^7 E A ls- / / ( £ ) Im [G (*, £ - en) G (k, -e + e kl ) +G(k,e + e kl ) G (k, -e - e kl )] de, (36) 



and 



A * = ~Tn S e * k ' hA Is- / / ( £ ) Im [G ( fc > e-e kl )G (k, -e + e kl ) + G(k,e + e kl ) G (k, -e - e kl )} de, (37) 



where 



G(k,e) = 7 -^- k , E k = ] /st Q + (Ali)\ (38) 



T 1 T 1 

_ _ £ k + £ k _ _ £ fc ~ £ fc 

£ /co — — ~ — , em. — — ^ — , i-jyj 



and Alg is given by relations (fT8|) and (fT9|) . The Fock's parameter appearing in Eqs. (Tl9|) is the Fourier transformation 
of I a = (cf a Cj a ). Based on the Zubarev's relation it can be calculated as 

^ = "a^ E (\ E elk ' h ) ^ / / ( £ ) Im (« c ^T : 4 T »e - ((d fcT ; c_ fcT » £ ) de. (40) 
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The above equations together with the equations for electron concentration and magnetization 



n=1 + ^E^ / /(£)Im((( CfcT ; C +)) £ + << Cfci ;c+)) £ -(( C i fcT ; C _ fcT )) e -(( C + fei ; C _ fei )) e )d e , (41) 
k J 



k " 



(42) 



constitute the set of self-consistent equations for parameters of ferromagnetic and superconducting states. 
The parameters Aq and A^, described by Eqs. ([36]) and (|37|) . can be calculated by the help of identity 



P 



1 



£-£k 



iirS (e - e k ) , 



(43) 



where P is the principal value of the integral and the symbol + stands for the infinitesimal positive value. As a result 
we obtain 



A = TT7 



1 \ - A n 

— ^£ [ tanh 
E k 



E k + £fei , , , Ek — Efcx 
tann 



2k B T 



2k B T 



(44) 



A 



J_ \ " ik-h 
AN ^ 



Alt / E k + e kl E k - eh 

kS 1 tanh „ + tanh k kl 



E k 



2k B T 



2k B T 



(45) 



In a similar way we obtain from the Eqs. (|4ip. (|4"2")l . and (|4H)l the following relations 



2N ^ 



£ ko I . , E k + e k i E k — £ki 

tanh — — — — h tanh 



E k 



2k B T 



2k B T 



(46) 



m 



E| tanh 



2N 



E k + £fei E k — e k i 

tann 



2fc s T 



2fc s T 



(47) 



J CT = - J- Y ( 1 V e* k ^ ( tanh ^*±££* _ Et^lMZ tanh gfc - g£fcl ) . (48) 

2N^\z^ j\ 2E k 2k B T 2E k 2k B T ) V 1 

Inserting Eqs. (|44f and (|45f into Eqs. (| 18|) and (|19p . and using the above relations for carrier concentration, magne- 
tization, and Fock's parameter, we can analyze the critical temperature versus carrier concentration for the singlet 
superconductors. The results are shown in Fig. 1. For simplicity it was assumed that 



la = n a (1 - n a ) , (49) 

instead of Eq. (|4"8"| . This result is strict for the constant density of states (DOS) at zero temperature. 

The critical temperatures shown in Fig. 1 depend on all interactions appearing in the Hamiltonian (|TJ) . In these 
and future calculations we use the dispersion relation for the 2 dimensional simple cubic lattice. 

As can be seen from this figure, the superconductivity created by kinetic interactions is not entirely pushed away 
by ferromagnetic ordering. In the interval of carrier concentrations around and below 1.5 superconductivity and 
ferromagnetism coexists together. The main driving force for the singlet superconductivity is the hopping interaction 
At. This interaction creates superconductivity by contributing to the negative pairing potential and also by decreasing 
the width of both spin sub-bands. The decrease of the majority spin sub-band width causes a gain in kinetic energy, 
which in turn creates a ferromagnetic state [22| . Such a ferromagnetic state, created by bandwidth changes, can 
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FIG. 1: Dependence of critical superconducting temperature for the singlet superconductivity (solid line) and Curie temperature 
(dashed line) on carrier concentration n. Dotted line shows the critical SC temperature vs. carrier concentration without 
ferromagnetism. Calculations were performed for J = 0.0875D, U = 0.75D, and D = 0.75 eV. The hopping inhibiting factors 
are: Si = 0.5 and S 2 = 0.05. 



coexist with the superconductivity |18j. This case is opposite to the effect of the band shift (i.e. Stoner model), which 
creates ferromagnetism but destroys superconductivity. 



4. Coexistence of ferromagnetism and triplet equal spin pairing (ESP) 
superconductivity 



The experimental evidence shows 36|, |37| ] that only the triplet parallel spins superconductivity can coexist with 
ferromagnetism within the same band of electrons. To explain this effect we assume in Eq. pip that the opposite 
spins parameters are zero (A^ = and A^ = 0) and we obtain the following relation 



£ — Si 



2(Ai)* e + el k 
V 2(Aj:)* 



2AJ 


He 




2Ai 



G(fc,e) = i, 



(50) 



where the energy gap A£ is given by Eq. (f2"5)) . 

Solving Eqs. (|50|) and (|25[) in a similar way to the previous case, we obtain the following equations for the supercon- 
ductivity 



A ^ = ~Jj E [ f (e) \G a (k, e) G" (k, -s)} de, 



(51) 



where 



G°{k,e) = {e-EZ)-\ E% 



+ (2A^) 2 . 



(52) 



In the H-F approximation, the superconductivity equation takes on the following form 



A£ = ly e * k - h ^itanh E? 
N ^ E? 



k 



2k B T' 



(53) 
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FIG. 2: Dependence of the superconducting critical temperature Tsc an d ferromagnetic Curie temperature Tc on carrier 
concentration n in the case of ESP for J = 0.0875D, Si = & = 0.4, U = 0, and D = 0.5 eV. Dashed line - Curie temperature; 
solid line - superconducting critical temperature without ferromagnetism; dot-dashed line - Tsc for Ai phase; dotted line - 
Tsc for A2 phase. 



From expressions (|4"Tj) . (|42p. and (|4T)]) the equations for carrier concentration, magnetization, and Fock's parameter 
can be calculated as 



n =i--y 

N 4^ 



2EI 



tanh 



■ tanh 



E, 



2k B T 



(54) 



— — tanh ^ fc 

^VV 2 4 2keT 



E, 



-^r tanh , " 
2E\ 2k BT/ 



(55) 



/ CT = — Y fiy e * h Vl-— tanh-^-V (56) 
2iV^-^^ J \ El 2k B T) K > 

Equations (|55|) - ([5rj)) are a set of self-consistent complex equations describing the ESP superconductivity coexisting 
with ferromagnetism. From this set we can find the dependence of the critical temperatures on carrier concentration. 
The results are shown in Fig. 2 for the Ai phase (where only parallel to magnetization pairing occurs) and for the A2 
phase (where antiparallel to magnetization pairing exists). 

For the superconducting effect without ferromagnetism the dependence of critical superconducting temperature on 
carrier concentration is determined by two factors: pairing potential c/4 given by Eq. (|26p and the change in kinetic 
energy described by the bandwidth modification factor W (n). Pairing potential is symmetric with respect to 
n = 1, where it reaches its maximum value. For the parameters J = 0.0875D, 51=5*2 = 0.4, and D = 0.5 eV we 
obtain dV jdn < at n < 1.6 and db a jdn > at n > 1.6. The decrease of b a factor causes an increase of the critical 
temperature. The net result of both these effects gives the maximum of Tsc{ n ) around n ks 1.25. When the inter-site 
correlations are not taken into account, V (n) = const, then the dependence Tsc(n) is determined solely by pairing 
potential. As a result the maximum of the critical temperature is localized at n = 1 (see e.g. [25j]). 

At concentrations with nonzero magnetization: n < 1.23, we have different Tsc f° r Ai and A2 phases. This 
difference in Tsc{n) is determined mainly by db a /dm. For A2 phase the increase in magnetization causes a decrease 
of the bandwidth modification factor: db^/dm < (see Fig. 3), that in effect increases the critical temperature. For 
Ai phase we have: dW /dm > 0, and the decrease in Tsc- The second factor favoring A2 phase at n > 1 is the shift 
of Fermi level for the minority carriers towards the center of the band, where there are higher critical temperatures 
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FIG. 3: Dependence of the bandwidth modification factor b" on magnetization m for n = 1.23, J = 0.0875D, and Si = Sa = 0.4. 
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FIG. 4: Dependence of the superconducting critical temperature Tsc for Ai and A2 phase on the Coulomb interaction U in 
the case of ESP for n = 1.23, J = 0.0875D, Si = S 2 = 0.4, and D = 0.5 eV. 



The difference between critical temperatures for Ai and A2 phases increases with the magnetization. This effect is 
visible in Fig. 4. Increasing Coulomb interaction U creates a growing exchange field in the H-F approximation, which 
induces increasing magnetization. This behavior was already reported by us [13]. Spalek and co-workers [H} have 
shown a similar effect of change in the superconducting gaps for Ai and A2 phases under an applied field. 



5. Influence of external pressure on superconductivity and ferromag- 
netism 

In the ferromagnetic ZrZri2 compound the ambient pressure level affects the Curie temperature. This compound 
has the quasi-linear dependence of both magnetic moment and Curie temperature on pressure. The experimental 
data are reported in Ref. [3§|. 

In addition to the quasi- linear dependence there is a rapid drop of the magnetic moment at pressure p c = 16.5 kbar, 
which seems to be the first-order phase transition. This transition may be caused by coupling between the long 
range itinerant magnetization modes and the weak particle-hole excitations. The coupling causes appearance of the 
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non-analytical terms in free-energy near the phase transition point [40| ■ 

In the ZrZn2 compound the superconducting phase coexists with ferromagnetism. For small pressures the super- 
conducting critical temperature is about 100 times smaller than the Curie temperature. 

To explain the pressure dependence of ferromagnetism and superconductivity we assume that the inter-site constants 
t in the Hamiltonian ([1]) depend on the external pressure. As a result, the kinetic interactions At and t ex will depend 
on pressure. They depend on t(p) through Eq. ^fy and on the hopping inhibiting factors S± and £2, which are pressure 
dependent. The on-site Coulomb repulsion U is assumed to be pressure independent. 

To find the pressure dependence of the hopping integral we need the results for the effective mass m* dependence 



on pressure [41 1 



m* (0) dp 
hence 



1 dm*(p) = _ 0017±0004kbar -i ) (57) 



m * (p) - 1.03(1 -Ap), (58) 



m 



where mo is the free electron mass and A « 0.017 ± 0.004 kbar -1 . 

Comparing the dispersion relation in the tight binding approximation at small k with the expression for nearly free 
electrons: = h 2 k 2 /2m*, one can obtain 



2m*a 2 2rriQa 2 m* 

Assuming for ZrZn2 that lattice constant [36| a = 7.393 A and inserting Eq. |55j) to Eq. (p)9")) one obtains the relation 



41] 



, , t(Okbar) 

m = \^. (60) 



We assume the following pressure relations for the hopping inhibiting factors S\ and S 2 



Si (p) = S 2 (p) = 1 
51(0 kbar) £2 (0 kbar) I -Bp 1 ' 

with parameter B smaller than parameter A. 

Pressure dependence of: t, t\, t 2 , and of the kinetic interactions: At and t ex , is presented in Figs. 5(a) and 5(b), 
for parameters A — 0.013 kbar" 1 and B = 0.003 kbar" 1 . 

Figures 6 and 7 show the pressure dependence of the Curie temperature. Figure 6 shows the dependence of Tc(p) 
for different carrier concentrations. Hopping integral t(p) is expressed by Eq. (|60[) and inhibiting factors S\ and S 2 are 
assumed to be pressure independent. The results show a decrease of the Curie temperature with increasing pressure. 
The Curie temperature always drops to zero with pressure. The exception is at the half-filling concentration. The 
explanation is following: pressure increase causes an increase of the hopping integral which in turn causes an increase 
of the critical (minimal) exchange interaction necessary for ferromagnetism while the existing exchange interaction 
(U) remains constant. In Fig. 7 we show the dependence of Tc(p) for different values of parameter B in Eq. (fCTj) . At 
positive values of B the decrease of Tq with pressure is faster than at B = 0; at B < we may have even an increase 
of Tc with pressure. 

Papers [la, EH show that ZrZn2 has a triplet parallel spin SC. Therefore the calculations will be performed 
for the coexistence of F with triplet SC. They will be based on Eqs. ([53 |) -(f55" f which give the triplet equal spin 
superconductivity gap, carrier concentration, and magnetization, respectively. The Fock's parameter will be expressed 
by the simplified Eq. (|4"9"|) . At critical temperature we assume A£ = in these equations. In Fig. 8 we present the 
dependence of Tsc(p) in the A2 phase at different values of parameter B. As shown in the previous Section, the 
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FIG. 5: Dependence of (a) hopping integrals and (b) kinetic interactions on pressure. The following values were used: A 
0.013 kbar -1 , B = 0.003 kbar -1 , Si (Okbar) = 0.35, S 2 (Okbar) = 0.25, and t (Okbar) = 0.125eV. 



critical temperature in the A2 phase is higher than in the paramagnetic phase. For positive values of B the critical 
temperature decreases with pressure, while for B < it increases with pressure. 

Figure 9 shows the dependence of Curie (Tc) and superconducting (Tsc) critical temperatures on the pressure. 
We used the parameters: A — 0.013 kbar -1 and B = 0.003 kbar" 1 , for which the hopping integrals and kinetic 
interactions increase with growing pressure [see Figs. 5(a) and 5(b)]. On the other hand, the kinetic interactions 
expressed in the units of bandwidth: t ex /t and At/t, decrease with growing pressure [see Fig. 5(b)]. In effect, both 
critical temperatures for superconductivity and ferromagnetism (Tsc an d Tc) decrease with pressure. This result is 
compatible with what was observed in ZrZna. 



6. Conclusions 

We have studied the possibility of coexistence between ferromagnetism and superconductivity within the framework 
of the extended Hubbard model. The main driving forces for these phenomena are the kinetic interactions: hopping 
At and exchange hopping t ex interaction. We have analyzed these interactions using the mean-field approximation. 

It was found that the hopping interactions influence ferromagnetism in two ways: 

- by increasing the Stoner (exchange) field, which shifts the spin sub-bands with respect to each other. This effect 
depends on carrier concentration. There is no half-filled band symmetry (n = 1). The effect is strongest for the 
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FIG. 6: Dependence of Curie temperature on pressure for different values of carrier concentrations. The parameters are: 
A = 0.017kbar _1 , Si = 0.35, S 2 = 0.25, and t (Okbar) = 0.125eV. 
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FIG. 7: Dependence of Curie temperature on pressure for different values of parameter B in Eq. (|61[) . The other parameters 
are: A = 0.013 kbar" 1 , Si = 0.35, Si = 0.25, and t (Okbar) = 0.125eV. 



smallest carrier concentrations and its strength decreases with increasing carrier concentrations; 

- by modifying the bandwidth factor b a , which depends on magnetization and carrier concentration. Decrease of 
W increases DOS, which favors ferromagnetism. This is a mechanism for raising the ferromagnetism from a gain in 
kinetic energy rather than from a decrease in potential energy. The effect is in favor of ferromagnetism for carrier 
concentrations n > 1. 

The kinetic interactions are also driving superconductivity. Singlet s-wave superconductivity is driven by the 
hopping interaction At, while the triplet superconductivity is driven by the exchange hopping interaction t ex . 

Analyzing the coexistence of ferromagnetism and the s-wave superconductivity we find that even weak ferromag- 
netism, if generated by the band shift, destroys the superconductivity. On the other hand the ferromagnetism created 
by a change of bandwidth can coexist with the singlet superconductivity. 

In the case of triplet superconductivity the ferromagnetism creates different critical SC temperatures for the Ai 
phase (where the pair's spin is parallel to magnetization) and for the A 2 phase (where the pair's spin is antiparallel 
to magnetization). With increasing magnetic moment the difference between both critical temperatures grows. This 
difference is caused by different bandwidths and different location of the Fermi level in spin sub-bands. 

The influence of pressure on magnetism and superconductivity was analyzed based on the assumption that we 
have: t(p), Si(p), and S^ip)- An increase of kinetic interactions induced by pressure causes a decrease of the Curie 
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FIG. 8: Dependence of superconducting critical temperature on pressure for different values of parameter B. Solid lines - Tsc(p) 
at nonzero magnetic moment, dashed lines - Tsc(p) in the paramagnetic state. The other parameters are: A = 0.013 kbar -1 , 
Si = 0.35, S 2 = 0.25, and t (Okbar) = 0.125 eV. 
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FIG. 9: Dependence of ESP superconducting critical temperature lOTgc and Curie temperature Tc on pressure. The parameters 
used are: n = 1.015, A = 0.013 kbar" 1 , B = 0.003 kbar -1 , Si (Okbar) = 0.35, S 2 (Okbar) = 0.25, and t (Okbar) = 0.125eV. 



temperature and also a change of critical temperature for triplet superconductivity. 
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